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INTRODUCTION 

For  a long  time  it  has  been  known  that  large  under- 
ground nuclear  explosions  have  small  Ms  values  compared 
to  earthquakes  having  the  same  mp  values  as  those  of  the 
explosions.  This  conclusion  has  been  found  also  to  be 
generally  true  for  Intermediate  and  amaller  sized  explo- 
sions. That  is,  when  a plot  of  the  mp  values  as  abscissae 
and  Na  values  as  ordinates  is  made,  the  earthquakes  separate 
from  the  explosions,  with  the  latter  lying  below  and  to 
the  right  of  the  earthquake  points  However,  it  has  been 
found  that  some  earthquakes  plot  in  or  near  the  explosion 
population.  These  are  called  anomalous  earthquakes. 

The  primary  purpose  of  this  research  Is  to  search 
for  shallow-depth  Eurasian  earthquakes  that  have  anoma- 
lous mb: values.  Furthermore,  because  it  is  for  inter- 
mediate to  small  magnitude  events  that  It  Is  most  diffi- 
cult to  discriminate  between  earthquakes  and  explosions, 
events  were  selected  which  had  mb  values  in  the  range  of 
4 to  5.5. 

Other  aims  of  the  research  are  to  identify  geographi- 
cal regions  where  anomalous  earthquakes  occur,  to  deter- 
mine the  effect.  If  any,  of  focal  depth  and  focal  mechanism 
on  observed  mb:  Mg  values,  to  Investigate  differences  in 
the  P-wave  spectra  of  anomalous  and  non-an omaloua  earth- 
quakes, and  to  determine  Love-wave  Mg (Mg  L)  values  for 
anomalous  earthquakes.  ' 


METHODOLOGY 

Standard  techniques  were  used  for  determining  the  ah, 
and  Mg  values  of  the  earthquakes  and  explosions  studied. 
That  Is,  for  mb  evaluation  the  largest  amplitude  in  the 
first  3 cycles  of  the  P-wave  motion  on  the  vertical- 
component,  short-period  seismogram  was  used,  along  with 
the  Gutenberg-Richter  (Richter,  1958)  calibrating  func- 
tion. Average  nip  values,  along  with  their  standard 
deviations  and  the  number  of  stations  used,  are  reported, 
as  well  as  average  mb  values  and  the  number  of  stations 
as  given  by  the  National  Earthquake  Information  Center 
and  the  International  Seismologlcal  Centre . 

To  determine  Mg  the  largest  amplitude  of  the  verti- 
cal component  Rayleigh-wave  motion  in  the  period  range 
of  17  to  23  seconds  was  used,  along  with  the  formulas 
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Us  • 3*30  + 1.66  log  A + log  A/P 


and 


for  25°  < A * 140° 


M3  ■ 4.16  + 1.07  log  A + log  A/P 


for  10°  £ a * 25° 


where  A is  the  eplcentral  distance  In  degrees,  A Is  the 
isaxlaua  ground  notion  In  microns,  and  T is  the  period 
In  seconds.  We  examined  the  use  of  Me  foramlas  derived 
for  waves  of  periods  substantially  less  than  20  seoonds 
as  reoorded  at  snail  to  regional  distances,  and  found  in 
general  that  they  were  not  as  satisfactory  as  the  seoond 
of  the  fornulas  given  above  (Huttll  and  Kin,  1975).  Mean 
Ms  values,  along  with  their  standard  deviations  and  nunber 
of  stations  used,  are  presented,  in  general  the  National 
Earthquake  Information  Center  does  not  determine  Ms  values 
for  snail  to  lnternedlate  magnitude  events,  and  the  Inter- 
national Selsmologloal  Centre  never  determines  Ms  values. 


The  seismograph  stations  whose  data  were  used  to 
determine  n^  and  Ms  values  were  the  World-Wide  Standard 
Selsnograph  Network  (WWSSN)  stations  ALQ,  AQU,  BLA,  BUL, 

££!'  25'  12k*  K0D*  K0M'  HIL>  lps,  mat, 

KBS,  NDI,  NBA,  00D,  QUE,  SEO,  SKI  and  SF&  and  the  Very 
Long  Period  Experiment  (VLPE)  stations  FBK,  ALQ,  TLO, 
CW,  CTA  and  EIL. 


Pocal  depths  were  determined  for  a selected  nunber 
of  earthquakes,  principally  by  means  of  the  depth  phases 
pP  and  sS.  When  depths  of  earthquakes  in  an  aftershook 
sequence  were  compared,  we  also  used  the  period  of  the 
minimum  in  the  Rayleigh-wave  spectrum . The  latter 
method  is  principally  of  value  In  establishing  that  two 
or  more  earthquakes  of  the  same  region  have  the  same 
depth,  rather  than  in  making  an  absolute  determination 
of  that  depth. 


DATA 

Much  of  the  basic  m^Mg  data  have  already  been  pre- 
sented in  Scientific  Report  No.  1 (Nuttll  and  Kim,  1974) 
and  Scientific  Report  No.  3 (Nuttll,  Kim  and  Wen,  1975). 
In  Report  No.  l data  for  96  earthquakes  and  10  under- 
ground explosions  were  given,  and  in  Report  No.  3 the 
®bMs  values  for  155  earthquakes  and  4 explosions  were 
presented.  Inasmuch  as  these  reports  are  readily  avail- 
able, a relisting  of  these  data  will  not  be  given  here. 
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Table  1 presents  the  hypocentral  coordinates  of  37 
earthquakes  and  9 explosions  whose  mb:  Mq  values  have  not 
been  Included  in  any  previous  report.  Table  2 gives  the 
aq>  and  Mg  values  of  the  events  described  In  Table  1. 

The  23  nuclear  explosions  took  place  at  the  testing 
sites  In  Novaya  Zemlya,  eastern  Kazakhstan  and  western 
Kazakhstan  and  at  scattered  sites  In  western  Russia. 

Figure  1 Is  a plot  of  the  mb: Mg  values  of  these  explosions. 
The  dashed  line  In  the  figure  Is  the  curve  obtained  by 
Evernden  et  al  (1971)  from  Nevada  Test  Site  data  when  Mg 
was  determined  from  20-second  period  Rayleigh  waves.  The 
two  solid-line  ourves  mark  the  envelopes  or  outer  bounds 
of  the  Eurasian  data.  From  the  upper  of  these  two  ourves 
we  can  see  that  all  the  explosion  data  satisfy  the  rela- 
tion 

®b  - Mg  >_  1.2  for  mp  > 4.6. 

The  single  point  n»b  - 4.4,  Ms  - 3-3  of  Figure  1 suggests 
that  for  mb  < 4.6  the  difference  mv.  - Ms  may  become 
smaller  than  1.2,  with  the  difference  decreasing  as  bh 
decreases.  However,  much  more  data  in  the  small  magni- 
tude range  than  we  have  are  required  to  Justify  such  a 
conclusion. 

Table  3 contains  a list  of  33  earthquakes,  out  of 
the  278  studied,  which  had  mb  minus  Mg  values  of  1.2  or 
greater,  and  which  thus  would  overlap  the  explosion  popu- 
lation on  an  mb: Mg  diagram.  The  mb:Ms  values  of  these 
earthquakes  are  plotted  in  Figure  2,  along  with  the  ex- 
plosion bounds  as  determined  from  Figure  1.  From  Figure 
2 It  oan  be  seen  that  the  anomalous  earthquakes  listed  In 
Table  3 have  mb:  Ms  values  that  spread  over  the  range  of 
explosion  values. 

Considering  next  the  non-anomalous  earthquakes,  our 
data  Indicate  that  shallow  oceanic  margin  earthquakes 
(Aleutians,  Kamchatka,  Japan,  Taiwan,  Philippines,  New 
Hebrides,  New  Britain)  all  have  n»b:  Mg  values  to  the  left 
and  above  the  upper-bound  curve  of  the  explosion  popula- 
tion. Thus  they  present  no  problem  in  being  distinguished 
from  explosions . Figure  3 1®  • plot  of  aq>:  Mg  values  of 
oceanic  margin  earthquakes. 

All  of  the  anomalous  earthquakes  that  are  Hated  In 
Table  3 and  plotted  in  Figure  2 are  Intraplate  earthquakes. 


J 





*p 


S95 


V'/ 

\ 


1 

('  \ 


Not  all  Intraplate  earthquakes , however,  are  anomalous. 
Thus,  for  example,  of  the  128  Intraplate  earthquakes  that 
we  studied  for  the  year  1972,  only  23  had  anomalous  mp:  Mg 
values . 


For  the  1972  earthquakes,  surface-wave  magnitudes  were 
calculated  using  Love-wave  amplitudes  as  well  as  the  ampli- 
tudes of  the  vertical-component  Rayleigh  waves.  Of  the 
122  earthquakes  for  which  Love-wave  M3  values  oould  be  de- 
termined, only  19  had  |Ms,r  - «s,l!  > 0-3-  Thus,  in  most 
oases,  the  Mg  » and  Mg  l values  were  essentially  the 
same.  Of  the  19  which' showed  a significant  difference, 
only  6 had  Mg  r > Ng^  R . Thus  there  is  little  evidence  of 
small  Mg  d values  compared  to  Ms  r values,  as  might  be 
expected  to  occur  if  the  focal  a£pth  were  such  that  there 
would  be  a minimum  in  the  Rayleigh-wave  spectrum  at  a 
period  near  20  seconds.  From  this  we  conclude  that  focal 
depth  (unless  it  is  very  large,  l.e.  greater  than  the 
crustal  thickness)  cannot  generally  be  called  upon  to  ex- 
plain small  Ms  values  of  earthquakes,  and  thus  anomalous 
mb: Mg  values. 


Evernden  (1976)  discussed  the  use  of  Mg  L 

58l6r 


— — values  to 

discriminate  between  earthquakes  and  explosl6ns.  At  plaoes 
where  the  release  of  regional  tectonio  strain  by  the  ex- 
plosion is  low,  the  resulting  Mg  r values  are  small  com- 
pared to  the  Mg  o values.  Of  the  23  explosions  considered 
in  this  report.  Love  waves  of  20-second  period  oould  be 
identified  for  only  two.  The  first  la  the  East  Kazakhstan 
explosion  of  16  August  1972,  for  which  a very  small  ampli- 
tude Love  wave  was  tentatively  identified  at  a single 


is  3.7  and 

aya  Zemlya 
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VLPE  station  (KON).  For  this' event  Mo  R ii 
is  3*0.  The  second  event  occurred  in  Nova] 

28  August  1972.  It  was  a large  explosion,  and  20-second 
period  Love  waves  were  identified  at  18  stations.  For  it 
MS  R was  5.0  and  Ms  L was  4.6. 


MA INSHOCK- AFTERS HOCK  SEQUENCES 


One  of  the  most  Important  findings  of  this  study  was 
that  some  aftershocks  of  a sequence  can  have  anomalous 
mb:  Mg  values,  whereas  the  malnshock  and  other  aftershocks 
have  non-anomalous  values.  This  finding  is  particularly 
Important  beoauae  it  enables  us  to  obtain  a better  under- 
standing of  the  cause  of  anomalous  earthquakes,  by  elimina- 
ting as  an  explanation  certain  characteristics  that  are 
often  proposed  as  the  causes  of  anomalous  mb:Ms  values. 
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Among  these  are  the  source  and  receiver  crust  and  the 
transmission  path,  which  are  identical  for  all  events 
of  a sequence  for  a given  seismograph  station,  and  the 
focal  depth  and  focal  mechanism,  which  are  shown  to  be 
the  same  for  selected  events  of  selected  sequences. 

A more  complete  discussion  of  the  data  and  con- 
clusions is  given  in  Appendix  1,  which  is  a manuscript 
of  a paper  submitted  for  publication.  In  this  section 
of  the  report  we  merely  present  some  of  the  data  and 
summarize  the  conclusions. 


Nalnshock-aftershook  sequences  in  Szechwan  province 
of  China,  Tibet,  the  Iraq- Iran  border  region  and  the  New 
Hebrides  islands  were  investigated.  Table  4 lists  the 
hypocentral  coordinates  and  the  mt>  and  Ms  values  of  the 
earthquakes  studied.  For  the  events  which  were  of  suffi- 
ciently large  magnitude,  focal  mechanisms  and  P and  Ray- 
lelgh-wave  spectra  were  determined.  The  focal  mechanism 
solutions  verified  that  the  earthquakes  of  a given  sequence 
had  nearly  identical  mechanisms,  and  the  Rayleigh-wave 
spectra,  along  with  the  time  differences  pP-P  and  sS-S, 
established  that  the  events  of  a given  sequence  had  nearly 
identical  focal  depths. 

From  the  P-wave  spectra  it  was  found  that  for  earth- 
quakes of  the  same  Mg  or  seismic  moment.  No,  the  anoma- 
lous earthquakes  had  smaller  spectral  corner  periods  than 
the  non-anomalous  ones.  This  indicated  that  there  is  a 
relative  enrichment  of  the  short-period  part  of  the  spec- 
trum for  the  anomalous  earthquakes.  Observations  in  the 
time  domain  showed  that  for  earthquakes  of  the  same  wk 
the  duration  of  the  short-period  P-wave  motion  was  smaller 
(more  pulse-like)  for  anomalous  than  for  non-anomalous 
earthquakes,  even  though  the  maximum  amplitude  of  the  P- 
wave  motion  was  the  same. 

Although  we  have  no  conclusive  explanation  as  to  the 
cause  of  the  anomalous  aftershocks,  the  observations  in 
both  the  frequency  and  time  domains  suggest  they  are  re- 
lated to  the  time  history  of  the  fault  rupture  and  pos- 
sibly to  the  stress  drop. 

RELEVANCE  OF  THE  RESEARCH  TO  THE 
DISCRIMINATION  OF  UNDERGROUND  EXPLOSIONS 

One  point  of  immediate  significance  is  the  finding 
that  anomalous  earthquakes  do  not  occur  at  plate  margins, 
but  rather  only  within  the  Eurasian  plate.  (We  did  not 
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look  at  earthquakes  In  other  continental  areas . ) Because 
earthquakes  occurring  along  plate  margins  constitute  the 
great  aajorlty  of  earthquakes,  our  finding  suggests  that 
for  this  large  class  of  earthquakes  the  nb:MS  criterion 
will  readily  separate  explosions  froa  earthquakes. 

A second  point  of  significance  is  that  the  anoaalous 
earthquakes  of  Eurasia  are  not  restricted  to  any  single 
or  few  geographic  areas.  Table  3 demonstrates  that  they 
occur  anywhere  from  southern  Europe  to  eastern  Asia . 

Although  froa  theory  one  aight  expect  a minimum  In 
the  Raylelgh-wave  spectrum  at  20-seoond  period  for  a 
certain  range  of  focal  depths  within  the  crust,  and  thus 
anomalous  05:  Mg  values  for  these  earthquakes,  our  data 
suggest  that  for  actual  earthquakes  recorded  at  a suffi- 
cient number  of  stations  this  Is  not  an  Important  phenom- 
enon. Nor  does  the  character  of  the  fault  motion,  whether 
It  be  strike-slip  or  thrust  faulting,  lead  to  anomalous 
mb:  Mg  values  by  Itself  If  the  Rayleigh  waves  are  recorded 
over  a sufficient  range  of  azimuth. 

The  mainshock-aftershock  studies  Indicate  that  for 
some  cases,  at  least.  It  Is  not  focal  depth,  focal  mech- 
anism, source  crust,  crust-mantle  transmission  path  or 
receiver  crust  which  Is  the  cause  of  anomalous  mv,:  Mg 
values.  Whatever  is  the  explanation,  it  causes  the  P- 
wave  spectra  to  be  enriched  at  the  short  periods  for  the 
anomalous  earthquakes  by  shifting  the  spectral  corner 
period  (or  periods)  to  smaller  values.  Our  finding  that 
anomalous  m^zMs  earthquakes  fairly  often  occur  as  after- 
shocks of  Intraplate  earthquakes  suggests  that  It  might 
be  a difficult  task  to  distinguish  between  a small  magni- 
tude explosion  and  an  aftershock  of  an  Intraplate  earth- 
quake, if  the  explosion  were  deliberately  detonated  after 
the  malnshook. 
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FIGURE  CAPTIONS 


Figure  1.  values  of  Eurasian  underground  explo- 

sions . The  values  are  those  of  the  International 
Selsmologloal  Centre,  and  the  Ms  values  are  those 
determined  in  the  present  study.  The  solid-line 
curves  mark  the  upper  and  lower  bounds  of  the  explo- 
sion data.  The  dashed-line  curve  is  the  NTS  explosion 
curve,  where  Ms  1®  determined  from  20-aecond  period 
surface  waves,  as  given  by  Evemden  et  al  (1971). 

Figure  2.  mb?Mo  values  of  anomalous  Eurasian  earth- 
quakes, for  which  nib-Mg  >,1.2.  The  solid-line 
curves  are  the  upper  and  lower  bounds  of  the  ex- 
plosion data,  as  given  in  Figure  1. 

Figure  3.  m>»:Mq  values  of  Eurasian  oceanic  margin 
earthquakes  of  1972.  The  solid-line  curves  are 
the  upper  and  lower  bounds  of  the  explosion  data, 
as  given  in  Figure  1. 
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Figure  1.  a*,:  Mg  values  of  Eurasian  underground  explo- 

sions. The  values  are  those  of  the  International 
Selsaologlcal  Centre,  and  the  Hg  values  are  those 
determined  in  the  present  study.  The  solid-line 
ourves  mark  the  upper  and  lower  bounds  of  the  explo- 
sion data.  The  dashed-line  curve  Is  the  NTS  explosion 
curve,  where  M3  Is  determined  from  20-second  period 
surface  waves,  as  given  by  Evemden  et  al  (1971). 


APPENDIX  1. 

SPECTRAL  AND  MAGNITUDE  CHARACTERISTICS 
OF  ANOMALOUS  EURASIAN  EARTHQUAKES 

by 

So  Ou  Kim*  and  Otto  V.  Nuttll 


ABSTRACT 

A number  of  malnshook-aftershook  sequences  In  the  Eurasian 
Interior  oontain  some  aftershooks  whose  ■*,:)%  values  sire  olose 
to  those  of  underground  explosions.  This  paper  is  ooneerned 
with  a study  of  the  amplitude  speotra  of  the  P waves  and  Ray- 
leigh waves  for  earthquakes  of  those  nainshoek-aftershook  se- 
quences. It  is  found  that  for  any  given  sequence  studied  there 
is  little  if  any  variation  in  focal  depth  or  focal  neohanism. 
This  rules  out  variations  in  these  quantities  as  being  the  cause 
of  anomalous  Ms  values.  A study  of  the  P-wave  spectra  estab- 
lishes that  one  or  both  of  the  corner  periods  of  anomalous 
earthquakes  are  smaller  than  those  of  non-anomalous  earth- 
quakes of  the  same  moment.  Thus  the  cause  of  anomalous  mt>: Mg 
values  of  the  earthquakes  studied  is  a relative  enrichment  of 
the  short-period  portion  of  the  spectrum  of  the  anomalous 
events,  which  cannot  be  attributed  to  focal  depth  or  focal 
mechanism . 


* Present  address:  Seismograph  Service  Corporation,  Tulsa, 

Oklahoma . 
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INTRODUCTION 

A recent  study  of  m^  and  Mg  values  of  approxljsately  300 
intermediate -else  Eurasian  shallow-focus  earthquakes  revealed 
that  most  anomalous  earthquakes  (events  whose  mpiNg  values  are 
closer  to  those  of  explosions  than  to  those  of  ordinary  earth- 
quakes) occurred  in  the  continental  interior  rather  than  at  the 
oceanic  margins  (Huttli  and  Kim, 1975) • Additional  study  has 
established  that  the  majority  of  anomalous  events  occurred  as 
aftershocks,  or  in  some  oases  foreshocks,  of  non-anomalous  main- 
shocks.  Rot  all  foreshocks  or  aftershocks,  however,  were  found 
to  be  anomalous. 

The  validity  of  using  m^Ng  values  in  distinguishing  between 
earthquakes  and  explosions  has  been  explained  (Holnar  et  al, 

19^9;  Lleberman  and  Pomeroy,  1970)  in  terms  of  a number  of  fac- 
tors, namely  souroe-tlme  functions,  source  size,  properties  of 
the  travel  path,  source  mechanism  and  focal  depth.  Tsai  and  Akl 
(1971)  argued  against  the  first  two  faotors  as  being  the  primary 
oause  of  the  difference  in  m^:  Mg  values  for  earthquakes  and  ex- 
plosions. Rather  they  maintained  that  the  last  two  are  dominant 
after  corrections  due  to  differences  in  travel  path  are  taken 
into  aooount.  Everaden  (1975)  also  concluded  that  anomalous 
events  oan  be  explained  principally  by  differences  in  source 
meohanism  and  focal  depth. 

On  the  other  hand,  Kanauorl  and  Anderson  (1975)  demonstrated 
that  inter-plate  earthquakes  whloh  occur  along  or  parallel  to 
major  plate  boundaries  show  low  stress  drop  oonpared  to  intra- 
plate earthquakes  which  occur  within  the  plates.  Forsyth  (1975) 


20 


! 

concluded  that  anomalous  m^:  Mg  values  observed  for  a foreshock- 
mainshock-aftershock  sequence  in  the  Kirglz-Sinkiang  border  re- 
gion resulted  from  the  source  finiteness  or  the  rupture-time 
duration  rather  than  from  the  focal  depth  or  the  focal  mechan- 
ism. 

EARTHQUAKES  STUDIED 

In  this  paper  we  restrict  our  discussion  to  single  mainshock- 
aftershock  sequences  in  the  Szechwan  province  of  China,  in  Tibet, 
in  the  Iraq-Iran  border  region  and  in  the  New  Hebrides  Islands. 

Table  1 gives  the  latitude,  longitude,  origin  time  and  focal 
depth  as  determined  by  the  National  Earthquake  Information  Ser- 
vice for  the  earthquakes  whose  focal  mechanisms  and  P-wave  and 
Rayleigh-wave  spectra  were  determined  in  this  study.  Also  in- 
cluded are  the  m^  and  M3  values  as  determined  by  us. 

An  inspection  of  Table  1 will  show  that  the  hypocentral  co- 
ordinates of  the  earthquakes  of  a given  sequence  show  little 
variation.  The  time  differences,  pP  - P and  sS  - S,  when  they 
could  be  observed,  indicated  a nearly  constant  focal  depth  for 
all  the  events  of  a sequence.  Further  evidence  for  the  same 
focal  depth  is  given  by  the  periods  of  the  minimum  in  the  Ray- 
leigh-wave spectra,  which  for  a given  station  are  almost  constant 
for  all  the  events  of  a sequence.  Sample  Rayleigh-wave  spectra 
for  two  stations  for  the  Szechwan  and  New  Hebrides  earthquakes 
are  given  in  Figures  1 and  2,  respectively. 

For  all  four  earthquake  sequences  the  focal  mechanisms  of 
individual  earthquakes  of  a given  sequence  are  remarkably 
similar.  This  is  demonstrated  in  Figures  3a  and  3b,  in  which 
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solutions  based  on  the  sign  of  the  onset  of  the  lung-period  P 
motion  are  presented. 

By  seleoting  for  study  events  of  nainshoek-aftershook  se- 
quences which  have  the  saae  hypooentral  coordinates,  the  saae 
foeal  mechanism  and  the  saae  transaisaion  path,  we  eliainate 
aany  of  the  quantities  whioh  have  been  proposed  in  explanation 
of  anoaalous  a^Mg  values.  We  can  thus  restrict  our  attention 
to  quantities  such  as  source  finiteness,  rupture  velocity  end 
stress  drop. 

Although  we  shall  not  present  a physical  explanation  as 
to  why  soae  earthquakes  have  anoaalous  a^: Mg  values,  we  wish 
to  point  out  certain  ways  in  whioh  the  spectra  of  anoaalous 
earthquakes  differ  froa  those  of  non-anoaalous  ones,  be  shall 
leave  to  others  the  interpretation  of  these  spectral  differ- 
ences, and  aerely  note  that  differences  in  the  corner  periods 
and  slopes  of  body-wave  speotra  have  been  interpreted  in  teras 
of  s.ource  finiteness  (Hanks  and  byss,  1972;  Forsyth,  1975)  snd 
of  fractional  stress  drop  (Brune,  1970) . 

TIME-DOMAIN  OBSERVATIONS 

Direct  observation  of  the  selsaograas  led  us  originally  to 
a study  of  aa ins hock-afters hook  sequences  in  a search  for 
anoaalous  events.  That  is,  by  looking  at  the  seismograms  we 
could  see  that  events  which  produced  nearly  the  saae  P-wave 
amplitudes  at  a station  produced  noticeably  different  Rayleigh- 
wave  aaplitudes.  Figure  4 gives  soae  examples.  On  the  left- 
hand  side  are  the  selsaograas  for  the  Szechwan  earthquakes. 

The  upper  left  shows  the  P-wave  motion  at  two  stations,  POL 
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and  SHK,  for  the  three  events.  Maxima  amplitudes  In  the  first 
few  cycles  are  about  the  same  (the  onset  of  the  P aotlon  at  SHE 
for  event  no.  3 la  obsoured  by  the  Minute  narks;  its  peak  and 
trough  almost  touoh  the  traoe  above  and  below,  respectively). 

The  Rayleigh-wave  amplitudes  for  event  no.  3,  as  reoorded  at  ata- 

tlon  SHI,  are  noticeably  less  than  for  events  1 and  2. 

i 

The  p waves  and  Rayleigh  waves  for  the  two  Tibetan  events 
are  given  in  the  right-hand  side  of  Figure  4.  The  P-waves  are 
shown  for  stations  MUR,  STU  and  XTO  and  the  Rayleigh  waves  for 
station  KTO.  Although  the  duration  of  the  P waves  of  event  no.l 
is  greater  than  of  no.  2,  the  amplitudes  of  the  first  few  cycles 
of  P of  no.  1 are  almost  the  same  as  the  amplitude  of  no.  2, 
which  would  result  in  nearly  similar  mb  values  for  the  two  earth- 

I 

quakes.  The  surface-wave  amplitudes  of  event  no.  1,  however,  are 
much  greater  than  of  event  no.  2,  as  can  be  seen  from  the  lower 
right-hand  figure. 

Figures  5 and  6 show  as  examples  m^rMg  plots  for  the  Tibet 
and  Mew  Hebrides  earthquakes,  respectively.  The  points  labeled 
1 and  2 in  Figure  5 refer  to  the  Tibetan  earthquakes  no.  1 and  2 
in  Table  1 . The  other  points  refer  to  earthquakes  for  which 
®b:*S  values  could  be  determined  but  which  were  too  small  for 


spectral  analysis.  Curve  a.  in  Figure  3 is  a least-square  linear 
fit  to  m^Mg  data  from  64  non-anomalous  Eurasian  earthquake 
values,  and  curve  a a similar  fit  for  11  Eurasian  underground 
explosion  values.  The  data  sets  were  taken  from  Muttli  and  Kim 
(1975).  Curve  b is  a linear  fit  to  the  Tibetan  data  points,  in- 
dicated by  x's  in  the  figure.  Note  that  the  mains hook  (event 
no.  1)  as  well  as  two  of  the  other  Tibetan  events  would  be 
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considered  non-anomalous . Event  no.  2,  as  well  as  6 others, 
would  be  considered  suspicious  by  the  criterion.  The  two 

remaining  events  definitely  fall  in  the  explosion  population. 

Curves  a and  o of  Figure  6 are  the  same  as  described  for 
Figure  5.  If  event  no.  2 is  considered  to  be  the  mainshock, 
its  M3  value  is  even  larger  than  for  the  ordinary  Eurasian  earth- 
quake of  mb  = 5-9-  The  foreshock  (no.  l)  and  all  the  aftershocks 
also  scatter  around  curve  a.  Curve  t)  is  the  linear  fit  to  the 
eight  data  points  of  the  figure . Thus  none  of  the  New  Hebrides 
earthquakes  would  be  considered  to  have  anomalous  mb:  Mc  values. 

RAYLEIGH-WAVE  SPECTRA 

Rayleigh-wave  ground-motion  spectra  were  obtained  from 
analysis  of  the  long-period,  vertical-component  seismograms. 
Stations  were  selected  so  that  the  path  from  epicenter  to  sta- 
tion was  pure  continental  or  pure  oceanic . The  observed  spectra 
were  adjusted  for  the  effects  of  geometric  spreading  and  anelas- 
tic  attenuation,  to  reduce  them  as  if  the  seismograph  station 
were  1000  km  from  the  epicenter. 

Sample  spectra  are  shown  in  Figures  1 and  2.  Values  below 
10  sec  are  unreliable  due  to  both  small  instrument  magnification 
and  low  signal  amplitude.  The  spectra  can  be  used  for  two  pur- 
poses: l)  to  determine  whether  differences  in  log  spectral 

amplitude  at  20-sec  period  are  equal  to  differences  in  Ms  values, 
l.e.  to  test  the  equivalence  of  amplitude  differences  in  the 
frequency  domain  with  amplitude  differences  in  the  time  domain, 

2)  to  determine  the  seismic  moment  of  the  various  earthquakes, 
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by  comparing  speotral  levels  of  actual  earthquakes  with  those 
calculated  for  a unit-moment,  double -couple  source  at  the  appro- 
priate azimuth  for  a distance  of  1000  las. 

Table  2 presents  the  Mg  differences  and  the  differences  In 
the  log  spectral  amplitude  at  20-seo  period  for  the  11  events 
In  the  4 geographic  regions.  The  number  n refers  to  the  number 
of  stations  whose  spectral  amplitudes  were  used  to  obtain  an 
average  log1Q  A A.  In  general  the  differences  between  & Mg  and 
lo*10  ^ * are  of  the  order  of  0.1  to  0.2  Mg  units.  The 

exception  is  the  set  of  Sxechwan  earthquakes;  either  Mg  of  event 
no.  1 is  too  snail  or  its  20-seo  period  spectral  amplitude  Is 
too  large . Prom  studies  of  the  P-wave  spectra  we  believe  that 
the  Mg  value  of  this  event  is  about  0.3  units  too  large,  which 
corresponds  to  one  standsrd  devistlon. 

Table  3 compares  the  seismic  moment,  as  obtained  from  the 
Rayleigh-wave  spectrum,  with  the  body-wave  magnitude.  Prom  the 
table  it  can  be  seen  that  non-anomalous  earthquakes  of  a given 
body-wave  magnitude  have  larger  moments  than  anomalous  earthquakes 
of  the  same  magnitude,  sometimes  by  more  than  a power  of  ten. 


P-WAVE  SPECTRA 

P-wave  ground-motion  spectra  were  determined  from  the  LPZ 
and  3PZ  seismograms  for  stations  at  teleselsmlo  distances.  As 
they  were  not  normalised  to  any  particular  eploentral  distance, 
in  comparing  them  one  should  compare  the  spectra  of  an  Individual 
station  for  an  individual  mainshook-aftershook  sequence.  Plgures 
8,  9,  10  and  11  give  sample  spectra  for  a few  of  the  stations 
whose  data  were  utilized.  The  small  open  circles  refer  to 
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speotral  values  obtained  froa  analyela  of  the  short -period  sels- 
aograas, and  the  oontlnuoue  curve  to  apeetral  values  obtained 
from  the  long-period  selsaograas.  Because  of  low  slgnal-to- 
nolse  ratios  the  speetral  values  obtained  froa  the  long-period 
selsaograas  are  not  valid  for  periods  less  than  3 sec. 

The  spectra  can  be  approximated  by  three  stralght-llne  seg- 
ments, namely:  a flat  portion  extending  froa  an  Infinite  period 
to  a corner  period  called  Tq^,  a sloping  straight  line  extending 
froa  the  corner  period  T01  to  a second  corner  period  T^,  and  a 
third  straight  line,  of  slope  greater  than  the  second,  extend- 
ing froa  comer  period  Ti2  to  the  short  periods.  The  first  and 
second  of  these  curves.  Intersecting  at  period  T^,  are  super- 
imposed on  the  spectra  In  Figures  7 through  10.  Froa  the  figures 
It  can  be  noted  that  In  general  the  non-anoaalous  earthquakes 


have  larger  values  of  TQ1  than  the  anoaalous  ones. 

Figure  11  shows  simplified  P-wave  spectra  for  each  of  the 
four  earthquake  sequences,  as  dete rained  froa  the  selsaograas 
of  a single  station  for  an  Individual  sequence.  In  aost  oases 
the  part  of  the  speotrua  between  the  comer  periods  Tq^  and  T12 
has  a slope  of  1.3  to  2.0,  and  the  part  between  T12  •***  1-sec 
period  has  a slope  of  3 to  4.  Differences  In  speotral  level  at 
1-seo  period  are  In  general  In  good  agreeaent  with  a^  differ- 
ences as  obtained  froa  tlne-doaaln  aeasureaents  of  the  selsao- 
graas. Differences  In  speotral  level  at  the  long  periods  agree 
In  a qualitative.  If  not  always  an  exaot  quantitative,  way  with 
differences  In  selsalo  aoaent  and  Mg . 
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The  ourves  for  MUR-Tlbet  of  Figure  11  show  that  although 
the  two  earthquake*  have  the  saae  m^  (1-sec  spectral  amplitude), 
the  non-anoaalous  earthquake  (no.  1)  has  longer  corner  period* 
and  a larger  seismic  moment . The  curve*  for  TIK-Szechwan  of 
Figure  11  also  show  approximately  the  same  speotral  level*  at 
1 sec,  hut  lesser  corner  period*  and  long-period  amplitude*  for 
anomalous  events  no.  3 and  1.  Similarly,  anomalous  Xraq-Xran 
events  no.  2 and  3 have  about  the  saae  speotral  amplitudes  at 
1 sec,  but  lesser  comer  periods  than  the  non-anoaalous  event 
no.  1.  If  we  consider  the  seismlo  moment,  proportional  to  the 
long-period  speotral  amplitude,  to  be  a more  fundamental  property 
of  earthquakes  than  the  m^  value,  then  we  oan  conclude  that  for 
a given  seismlo  moment  an  anomalous  earthquake  has  smaller  cor- 
ner periods  and  a larger  1-seo  speotral  amplitude  than  a non- 
anoaalous  one. 

DISCUSSION  AMD  CONCLUSIONS 

We  have  shown  that  earthquakes  with  anomalous  g*>:  Mg  values 
occur  in  some  of  the  foreshocks  and  aftershocks  of  sequences 
in  Szechwan  Provlnoe,  Tibet  and  Iraq-Iran  earthquakes  of  the 
Sura sian  interior.  The  earthquakes  in  a given  sequenoe,  non- 
anoaalous  and  anomalous  alike,  have  similar  focal  depth  and  fo- 
cal mechanism.  Thus  neither  differences  in  these  fooal  param- 
eters nor  in  the  transmission  path  from  source  to  station  can 
be  oalled  upon  to  explain  the  anomalous  earthquakes. 

In  all,  we  considered  29  earthquakes  in  the  regions  men- 
tioned above  plus  8 in  the  New  Hebrides  region.  Although  there 


is  IOM  aablgulty  as  to  what  constitutes  sn  anoaalous  ■btNs 

▼slue,  ss  aany  ss  24  of  the  29  oan  be  considered  to  be  anoaa- 
lous  or  near-anoaalous  events.  All  of  the  6 Mew  Hebrides 
events,  on  the  other  hand,  are  non-anoaalous . 

We  found  that  a oossson  feature  of  the  snosialous  events 
was  that  the  corner  periods  of  their  P-wave  spectra  were  less 
than  those  of  non-anoaalous  events  which  had  the  saae  value. 
We  also  observed  In  the  tlae  domain,  that  for  earthquakes  of 
the  saae  a^the  duration  of  the  short-period  P-wave  notion  was 
less  (aore  pulse-llke)  for  the  anoaalous  than  for  the  jon- 
anoaalous  events. 

The  findings  of  this  study  have  tone  bearing  on  the  problea 
of  nuclear  test  detection.  In  addition  to  the  obvious  one  that 
soae  Intra-plate  foreshocks  and  aftershocks  possess  anoaalous 
BfetMg  values.  Because  these  anoaalous  values  oan  be  explained 
In  terns  of  shorter  corner  periods  and  thus  an  enrlohaent  of 
the  short -period  portion  of  the  p-wave  spectrua,  rather  than 
In  differences  In  the  long-period  portion  or  the  selsalo  aonent. 
It  follows  that  aethods  of  dlscrlnlnetlon  that  eoapare  the  long- 
period  portions  of  explosion  and  earthquake-generated  waves, 
such  as  relative  excitation  of  Rayleigh  and  Love  waves,  will 
be  unaffected  by  the  results  discussed  here.  Thus  the  aethods 
of  discrimination  which  utilise  exclusively  the  long-period 
portion  of  the  surface  and/or  body-wave  spectra  should  be  use- 
ful in  discriminating  between  explosions  and  those  relatively 
few  earthquakes  which  are  classified  as  anoaalous  by  the  conven- 
tional a^: Mg  analysis. 
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HYPOCENTRAL  COORDINATES  AND  MAGNITUDES  OF  EARTHQUAKES  STUDIED 
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TABLE  2 

COMPARISON  OF  AND  LOG  A_A  VALUES 

i Region  Earthquake  Pair  AMg  log  a A 


Szechwan  Province 

No. 

2 

- No. 

1 

0.6 

0.21 

(n 

= 

6) 

No. 

1 

- No. 

3 

0.1 

0.49 

(n 

= 

6) 

No . 

2 

- No. 

3 

0.7 

o 

o 

(n 

= 

6) 

Iraq-Iran  Border 

No. 

1 

- No . 

2 

0.5 

0.51 

(n 

= 

4) 

No. 

1 

- No. 

3 

0.8 

o 

00 

o 

(n 

= 

4) 

No. 

2 

- No. 

3 

0.4 

0.29 

(n 

= 

4) 

Tibet 

No. 

1 

- No . 

2 

1.20 

1.45 

(n 

= 

2) 

New  Hebrides  Islands 

No. 

2 

- No. 

1 

1.39 

1.26 

(n 

= 

4) 

No. 

3 

- No. 

1 

0.53 

0.37 

(n 

= 

4) 

No. 

2 

- No. 

3 

0.86 

0.99 

(n 

= 

4) 

I 
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TABLE  3 

COMPARISON  OF  mb  AND  M0  FOR  ANOMALOUS  AND  NON- 
ANOMALOUS  EARTHQUAKES 


Earthquake 

mb 

Mo 

Anomalous 

(dyne -cm) 

Non-anomalous 

Iraq-Iran  No 

. 2 

5.1 

5.5  x 1023 

Iraq-Iran  No 

. 1 

5-3 

1.9  x 

1024 

Iraq-Iran  No 

. 3 

5.3 

2.8  x 1023 

Szechwan  No. 

3 

5.3 

3-2  x 1023 

New  Hebrides 

No.  l 

5-3 

3 .6  x 

io24 

Szechwan  No. 

2 

5.4 

1 . 8 x 

1024 

Szechwan  No. 

1 

5.7 

1.1  x 1024 

New  Hebrides 

No.  3 

5-9 

8.5  x 

1024 

New  Hebrides 

No.  2 

5-9 

9.0  x 

1023 

Tibet  No.  1 

6.0 

3-3  x 

1023 

Tibet  No.  2 

6.0 

1.2  x 1024 
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Figure  Captions 


Fig.  1 


Fig.  2 


Fig.  3b 


Fig.  3b 
Pig.  4 


Rayleigh-wave  spectra  of  Szechwan  earthquakes  as  re- 
corded at  stations  NDI  and  QUE  and  normalized  to  a 
distance  of  1000  km. 

Rayleigh-wave  spectra  of  New  Hebrides  earthquakes  as 
recorded  at  stations  JCT  and  OGD  and  normalized  to  a 
distance  of  1000  km. 

Focal  mechanism  solutions  for  Szechwan  earthquakes 
(above)  and  Iraq-Iran  earthquakes  (below).  Compres- 
sions are  indicated  by  octagons,  dilatations  by  tri- 
angles and  small  amplitude  arrivals  by  X's. 

Focal  mechanism  solutions  for  the  Tibet  earthquakes 
(above)  and  the  New  Hebrides  earthquakes  (below), 
(above)  P-wave  motion  on  SPZ  records  for  the  Szechwan 
earthquakes  (left)  and  tne  Tibet  earthquakes  (right). 
For  the  Szechwan  earthquakes  the  station  on  the  left 
is  BUL,  at  an  epicentral  distance  of  87-3°  and  an  azi- 
muth of  19. 9° » and  the  station  on  the  right  is  SHK, 
at  an  epicentral  distance  of  25-9°  and  an  azimuth  of 
265. 2°.  For  the  Tibet  earthquakes  the  station  on  the 
top  is  NUR  ( A = 40.4°,  BAZ  = 94.1°),  in  the  middle  is 
STU  ( A = 50. 7°,  BAZ  = 73-9°)  and  on  the  bottom  is  KTG 
( A = 63.1°,  BAZ  = 0.0°) . 

(bottom)  Rayleigh-wave  motion  at  station  SHI  ( A = 44.3 
BAZ  = 77.8°)  for  the  Szechwan  earthquake  (left)  and  at 
station  KTG  ( A = 63.1°,  BAZ  = 60.6°)  for  the  Tibet 
eartnquake  (right). 


Fig.  5 Bb:Ns  diagram  for  the  Tibetan  earthquakes.  Noe.  1 and 

2 refer  to  the  earthquakes  listed  In  Table  1.  Curve  a 

Is  the  linear  least-square  fit  to  54  non-anomalous 

Eurasian  earthquakes.  Curve  b is  the  least-square  fit 

to  the  Tibetan  earthquake  data,  indicated  by  X's. 
e 

Curve  o is  the  least-square  fit  to  11  underground 
Eurasian  explosion  data. 

Fig.  6 n^iMg  diagram  for  the  New  Hebrides  earthquakes,  indi- 
cated by  X's.  Nos.  1,  2 and  3 refer  to  the  earthquakes 
listed  in  Table  1.  Curve  b Is  the  linear  least-square 
fit  to  the  New  Hebrides  data.  Curves  a and  o are  as 
in  Fig.  5. 

Fig.  7 P-wave  amplitude  spectra  at  stations  TIE  and  COL  for 
the  Szechwan  earthquakes.  The  continuous  curve  rep- 
resents the  spectra*  obtained  from  the  long-period, 
vert leal -component  seismogram  and  the  circles  the 
spectral  values  obtained  from  the  short-period,  ver- 
t leal-component  seismogram. 

Fig.  6 P-wave  amplitude  spectra  at  stations  STU,  NOR  and 
KTO  for  the  Tibetan  earthquakes. 

Fig.  9 P-wave  aaplltude  spectra  at  stations  AAE  and  SHL  for 
the  Iraq- I ran  earthquakes. 

Fig.  10  P-wave  aaplltude  speotra  at  stations  NAT  and  SHL  for 
the  New  Hebrides  earthquakes. 

Fig.  11  Simplified  P-wave  spectra  for  the  Sseohwan  earthquakes 
as  reoorded  at  station  TIK  (upper  left),  the  Iraq-Iran 
earthquakes  as  recorded  at  station  SHL  (upper  right). 
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the  Tibet  earthquakes  as  recorded  at  station  MUR 
(lomr  left)  and  the  New  Hebrides  earthquakes  as  re- 
oorded  at  station  SHL  (lower  right).  The  differences 


amplitude  at  1-seo  period  scale  as  differences  in 


*b»  And  the  differences  in  saplltude  at  20-sec  period 
*•  differences  in  Mg.  The  corner  period  T^  for 
Szechwan  earthquakes  no.  1 and  3 la  less  than  1 see. 
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